Benchmarking a semiclassical impurity solver for dynamical-mean-field theory: 
self-energies and magnetic transitions of the single-orbital Hubbard model 
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An investigation is presented of the utility of semiclassical approximations for solving the 
quantum-impurity problems arising in the dynamical-mean-field approach to the correlated-electron 
models. The method is based on performing a exact numerical integral over the zero-Matsubara- 
frequency component of the spin part of a continuous Hubbard- Stratonovich field, along with a 
spin-field-dependent steepest descents treatment of the charge part. We test this method by ap- 
plying it to one or two site approximations to the single band Hubbard model with different band 
structures, and comparing the results to quantum Monte-Carlo and simplified exact diagonalization 
calculations. The resulting electron self-energies, densities of states and magnetic transition temper- 
atures show reasonable agreement with the quantum Monte-Carlo simulation over wide parameter 
ranges, suggesting that the semiclassical method is useful for obtaining a reasonable picture of the 
physics in situations where other techniques are too expensive. 

PACS numbers: 71.10.Fd, 75.10.-b 



I. INTRODUCTION 

Correlated-electron materials, in which the interaction 
energy is comparable to or larger than the electron kinetic 
energy, are an important topic in materials science. In 
these compounds, standard band-theory is an inadequate 
representation of the physics. The discovery of high- Tc 
superconductivity in the oxide cupratesi led to greatly in- 
creased interest in correlated-electron compounds. Many 
materials have been studied, and many novel proper- 
ties have been discovered; including colossal magnetore- 
sistance, variety of spin, orbital, charge orderings, and 
unconventional superconductivityj^ Understanding the 
novel phenomena and determining the correct electronic 
phases in these materials are challenging tasks, and are 
indispensable for developing the electronic devices ex- 
ploiting the novel properties of the correlated-electron 
materials. However, the rapid increase of Hilbert space 
with system size limits the utility of exact-diagonalizaion 
(ED) methods, while the fermion "sign problem" renders 
Monte-Carlo (MC) approaches ineffective. 

Dynamical-mean-field theory (DMFT) is a promis- 
ing approach for treating correlation effects4i^ The 
method may be combined with realistic band-structure 
calculations, and has bee n applied to variety of ma- 
terials (see Refs. I6ll7ll8l9l for example). In DMFT, 
the momentum dependent electron self-energy is ap- 
proximated as a finite set of functions of frequency, 
and physical lattice problem is mapped into a one 
or several coupled impurity-Anderson models (lAM) 
with environment determined self-consistently. To 
solve the lAM, a variety of theoretical methods 
have been applied, including analytical approximations 
such as iterated-perturbation theoryji2iii*iS non-crossing 
approximatioufi^- projective self-consistent method,^'* 
and slave-bosoni^ for example, and numerical meth- 
ods including Quantum Monte-Carlo (QMC)fiSii2i exact 



diagonalization,^* numerical-renormalization group^SiSS 
and density-matrix-renormalization group^i*^ Most in- 
teresting phenomena can not be addressed by analyti- 
cal methods. Numerical methods, however, are compu- 
tationally very expensive. In order to combine DMFT 
with a realistic band-structure calculation and survey 
a wide range of parameters, it is desirable to develop 
computationally-cheap and reliable methods. Methods 
which can reproduce the correct electronic phases and 
distinguish the phases near in energy at low temperature 
are particularly needed. Many ideas have been proposed 
to reduce the computational costi^^i^'^i^^'^^i^'^i^* These 
fall mainly into two groups: variations of the method 
simplified by truncating the number of orbitalSf22i24 and 
perturbative methods involving expansions around a cer- 
tain solvable limit. ^j^*2SiSL2& The restricted method pro- 
vides a good approximation to the T = Mott physics 
of the Hubbard model^iS^ however because only a small 
number of states are used, it has difficulty dealing with 
the thermodynamics (as discussed later), and becomes 
impractical for multiorbital and multisite systems. Per- 
turbative methods are not reliable at intermediate cou- 
pling. Non-perturbative methods which can deal with 
the thermodynamics are required. 

In this paper, we investigate an alternative method 
to solve the lAM in DMFT formalism: a semiclassi- 
cal approximation based on the continuous Hubbard- 
Stratonovich transformation. ~ This method is compu- 
tationally very cheap (approximately two-orders of mag- 
nitude faster than QMC for single-impurity models, and 
with a better scaling with system size): it may be per- 
formed on a commercial PC. The semiclassical approxi- 
mation, in the form used here, was apparently first pro- 
posed by Hasegawa,'^*' who used it to study a "single- 
site spin fiuctuation theory" which may be viewed as 
an early, simplified version of dynamical-mean-field the- 
ory. Semiclassical methods were also used by Blawid 
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and Millie^"'' and by Pankov, Kotliar and Motome^-- to 
study models of electrons coupled to large-mass oscilla- 
tors. In this paper, we present an implementation in the 
context of dynamical-mean-field theory, and test its re- 
liability for a fully quantal model problem, namely, the 
Hubbard model on a variety of lattices by performing de- 
tailed comparisons between the semiclassical approxima- 
tion and QMC. We also present a brief comparison a suc- 
cessful "restricted method" : the two-site DMFTi2,3 The 
semiclassical approximation is found to be reasonable 
at high temperature and in the strong coupling regime, 
and gives good results for magnetic transitions. In the 
half-filled square lattice, the Neel temperature computed 
by the semiclassical approximation is very close to that 
found in QMC. In the metallic face-centered cubic (FCC) 
lattice, ferromagnetic Curie temperatures computed by 
the present method are within a factor of 2 compared 
with the QMC, and the competition between ferro- and 
antifcrromagnetism is correctly captured, T — phase 
boundaries are obtained within the error of ~ 20 % to 
QMC, The T = phase boundaries are in the same range 
as found in two-site DMFT, but this latter method gives 
very poor transition temperatures. We also use the semi- 
classical approximation to compute dynamical proper- 
ties, in particular, the electron self-energy and spectral 
functions, which are found to be in reasonable agreement 
with QMC results. 

The rest of this paper is organized as follows: Section II 
formulates the semiclassical approximation. In Sec. Ill, 
we use the Gaussian fluctuation approximation to esti- 
mate the limits of validity of the semiclassical approxima- 
tion. Sections IV and V compare numerical results for the 
semiclassical approximation to QMC, ED, and Hartree- 
Fock (HF) results, for the single-band Hubbard model 
on square lattice and face-centered cubic lattice, respec- 
tively. Section VI presents an application of the method 
to the two-impurity cluster dynamical-mean-field approx- 
imation for the half-filled Hubbard model. Comparison 
of nearest-neighbor spin correlation between the semi- 
classical approximation and QMC are given. In Sec. VII, 
we discuss the equilibration problem associated with the 
partitioned phase space in the semiclassical approxima- 
tion and QMC. Finally, Sec. VIII gives a summary and 
discussion. 



II. FORMULATION 

In this section, we derive the semiclassical approxi- 
mation from the continuous Hubbard-Stratonovich (HS) 
transformationiS We consider a single-orbital repulsive 
Hubbard model for simplicity: H = Hhand + Hu with 
Hband and Hu being single-particle dispersion term and 
interaction term, respectively. Generalizations to multi- 
orbital and long-ra nge interaction would be also possible: 
see e.g., Refs. Isals^lSSllSGl and WA 

We take the band-dispersion term to be Hband — 
J2kcr ^kcl^Cka with three choices of Sk- 



Two dimensional square lattice: 

e'j?^"'^'^ = —2t(cos kx + cos ky). (1) 
Three dimensional FCC lattice: 

e^^^^ — 4i(cos kx cos ky + cos ky cos k^ + cos k^ cos k^) 



-|-2i'(cos 2kx + cos 2ky + cos 2kz), 



(2) 



where t and t' are the NN hopping and the second NN 
hopping, respectively. 

Infinite- dimensional FCC lattice: Non-interacting den- 
sity of states is given by 



^FCCoo(^) ^ exp{-(l + V2£)/2} ^ 
U(l + V2e) 



(3) 



which diverges at the bottom of the band. ]S[P'^'^°° was 
introduced by Miiller-Hartmann4£ and studied in detail 
by Ulmkep^ The FCC lattices are of interest because 
these display a ferromagnetic phase in a wide range of 
doping. 

The interaction term is given by 



(4) 



with U >0. 

In the DMFT approximation, one first needs to com- 
pute the partition function of an A^-site impurity model 



with 



Z = J V[c\c]exp{-{So + S,nt)}, 
Jo 



(5) 



(6) 



where, ip = [ci, . . . , cat]' with Ci(c|) being a Grassmann 
number corresponding to the electron annihilation (cre- 
ation) operator at site or orbital i. is the N x N 
matrix Weiss field (inverse of the non-interacting Green 
function) which will be determined self-consistently, and 
/3 = l/T is the inverse temperature. Sint represents 
the interaction term specified by the model one con- 
siders. For the single-band Hubbard model, Sint — 
UJdTJ2in^^iT)n^liT). 

Next, one computes single-particle interacting Green 
function Go- by a functional derivative of InZ with re- 
spect to ao- as 



^In.^ 



(7) 



The electronic self-energy is obtained by inverting the 
Dyson equation as 



g: 



(8) 
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Finally, by connecting the impurity Green function G 
and the local part of the lattice Green function, the 
DMFT equation is closed. The self-consistency equation 
for the ij component of G is 



1.0 



27r 



/i 



(9) 

where, w„ is the fermionic Matsubara frequency, and tfc 
and Sfc are the Fourier transforms of the hopping and the 
self-energy matrices, respectively. (pij{k) is a form factor 
specified by the DMFT method choseni2& The chemical 
potential /i is fixed so that G gives the correct electron 
density n. The local Green function G is used to update 
the Weiss field a^'^'" as a^^"' = G^^^ -f S,,, and this pro- 
cess is repeated until the self-consistency a^r — a^^™ is 
obtained. The expensive computational task is evaluat- 
ing the functional integral in Eq. jSJl- 

A key point of the present approximation for the single- 
band Hubbard model is introducing the charge field as 
well as spin. Using the following identity: 



1 



(10) 



with p — n-^ + ni and m ~ — n^, we perform the HS 
transformation exp(— iSi„f) — J x] exp(— with 
the effective interaction 



1 /■'^ 



-2[/^4(T){(^(T)a, +zx(t)K(t) 



, (11) 



where, ip and x are the HS fields acting on spin and charge 
degrees of freedom, respectively, and az is the z compo- 
nent of the Pauli matrices. Now, one can perform Grass- 
mann integral over the fermionic field formally to obtain 
the partition function 



x] exp(— 5, 



(12) 



with the effective action 



1 

-Tr In 



/3 



dT{^^{T)+X^{T)} 



(13) 



Here, the trace includes spin as well as Matsubara fre- 
quency. It is noted that the coupling constant for the 
charge field x is imaginary. This originates from the dif- 
ferent sign of quadratic terms for spin and charge when 
decoupling the interaction term, [see Eq. 1(1171) ] 

Exact evaluation of the partition function Eq. \i2\ is 
impossible. The simplest approximation is to approxi- 
mate the integrals by the value computed using the static 
parts (p and x which extremize the action (solutions of the 
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FIG. 1: Example of the effective potential for ip, V{ip) given 
in Eq. in the square-lattice half-filled Hubbaed model 

with U/t = 6. Solid line: T/t = 0.4 (about 14 % above 
Tat ~ 0.35t); broken line: T/t = 0.3 (about 14 % below Tat). 
At both temperatures, V is far from a simple parabola, and 
regions far from the local minima have appreciable occupation 
probability. 



saddle-point equations dS{(f, x) jdCp — dS{ip, x)/dx = 0) . 
This is the Hartree-Fock (HE) approximation, which is 
known to give a poor estimates to transition tempera- 
tures and self-energies. One may correct the HE ap- 
proximation by including the Gaussian fiuctuations in 
which Sef / is expanded around the saddle point up to the 
quadratic order of the fluctuations 5ip{T) and 5x{t). In 
this case, 5(p{t) and 6x{t) decouple. However, Gaussian 
fluctuation theory is limited to the weak coupling regime. 
As an example. Fig. shows an effective potential V for 
(p, equivalent to TSeff, calculated (as explained below) 
by the semiclassical approximation. It is seen that the 
potential is highly non-parabolic, and the variation in V 
is on the scale of T for reasonable T. Thus, the Gaus- 
sian fluctuation theory is inapplicable. Next, one can 
consider taking tp and x only at Matsubara-frequency 
i^i ~ 0, i.e., static approximation, but evaluating the par- 
tition function Z as a two-dimensional integral over two 
static fields (2-field approximation). The partition func- 
tion is expressed as Zstatic ~ Jdipdxe'Kp{—Sstatic) with 



'^static — ' 



Trln 



(14) 



However, this approximation fails because the effective 
action Sstatic is complex and exp^—Sstatic) can not be 
regarded as a distribution function of cp and x, leading 
to poor convergence of integrals. With some effort, ap- 
parently converged integrals can be obtained, but the 
interacting Green function computed by using Eq. Q 
with Z replaced by Zstatic does not behave correctly; 
the imaginary part of the self-energy changes sign and 
the spectral-function sum-rule is strongly violated. As 
an example. Fig. [21 compares self-energies computed by 
evaluating the static average integrating over two static 




FIG. 2: Comparison of the electron self-energy computed by 
evaluating the static average integrating over two static fields 
ip and X with the action Eq. 11411 (2-field) and the semiclas- 
sical approximation, which is defined in this section, for the 
infinite-dimensional FCC lattice Eq. (|3J. Squares (circles) are 
real (imaginary) part of the self-energy, and filled and open 
symbols are obtained by 2-field approximation and the semi- 
classical approximation, respectively. 



fields ip and x with the action Eq. H14|l to the semiclassical 
approximation defined shortly. While we have obtained 
the converged solution along the Matsubara-axis, imag- 
inary part of the self-energy changes sign and causality 
is strongly violated, so that the analytic continuation to 
the real axis is impossible. 

In order to reduce the above mentioned problems and 
to take into account the fluctuation of both fields to some 
extent, we apply the semiclassical approximation follow- 
ing Hasegawat^S. In this method, we first solve the mean- 
field equation for the static charge field which extrem- 
izes 5 at a given value of ip. From dSeff /dx\^ = 0, we 
obtain the following equation (^^ = ix\ip), 



-UTTt- 



1 



(15) 



By solving Eq. p5|l . one obtains as a function of (p. 
In the single-impurity model, this equation has a unique 
solution, thus, Eq. IjlSfl can be solved easily, for example, 
via bisection. (For multiimpurity models, we have no 
proof that there is a unique solution, but difficulties have 
not arisen in the case we have studied.) Fluctuations of 
the charge field around this saddle point are expected to 
be less important than those of spin field, because of the 
unique solution and because charge fluctuations are not 
soft, and the difficulties associated with the oscillatory 
convergence of Eq. H14|) are avoided. 

The t/j-dependence of ^ is crucial. If it is neglected, i.e., 
if ^ is taken to be the value ^av which extremizes Seff 
after averaging over if, then one obtains a model equiva- 
lent to the Jahn- Teller phonon model studied for example 
in Ref. [s^ This latter model has different physics from 



FIG. 3: Comparison of the many-body density of states com- 
puted for the paramagnetic phase of the infinite-dimensional 
FCC lattice [Eq. Q] by the "full" semiclassical approxima- 
tion (solid line) to that computed using additional approxi- 
mation in which ^ is assumed to be independent of ip (broken 
line). The three-peak structure seen in the broken line is un- 
physical; the solid line is in good agreement with QMC (not 
shown). Note that Tc ~ 0.073 for these parameters; we have 
presented data at T = 0.05 ~ 3rc/4 (chosen to reveal the 
three-peak structure clearly), and have artificially suppressed 
ferromagnetism. 



the Hubbard model. In particular, the Jahn- Teller model 
has a spectral function with a characteristic three-peak 
structure quite unlike the spectral function of Hubbard- 
like models. An example of the density of states com- 
puted by the approximation, S^^p ^av is shown in Fig.|2| 
as the broken line. One can clearly observe a three-peak 
structure. The outer two peaks correspond to the oc- 
cupied (lower band) and unoccupied (upper band) state 
at the occupied distorted site (from potential minima at 
ip ~ ±?7), and the middle peak corresponds to the un- 
occupied undistorted state (from potential minimum at 
ip 0) in the phonon model. In the phonon model, the 
level separation between the lower occupied state and the 
middle unoccupied one has physical meaning as "polaron 
binding," which does not exist in the original Hubbard 
model. The spectral function computed from the "full" 
semiclassical method (see below) is shown as the solid 
line in Fig. 13 and as will be seen below is in much better 
agreement with QMC data. 

Now that the saddle point of the effective action is 
determined with respect to one of two variables, ^, an 
effective potential for remaining variable (p is written as 



(16) 



With this effective potential, the partition function is ap- 
proximated as 



dipeM-PVicp)}. 



(17) 
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It may be useful to mention the correspondence 
between the present theory and the single-site spin- 
fluctuation theory of Hasegawai^ In Ref. |^ Hasegawa 
introduced the same HS transformation using spin- and 
charge-fields, and applied saddle-point approximation 
against the charge-field. Thus, up to this point, two the- 
ories are equivalent, but instead of Eq. (|f 8|l and DMFT 
self-consistency equation, Hasegawa used an ad-hoc pro- 
cedure involving computation of averages of the magnetic 
moment and square of the magnetic moment. 

We emphasize, however, that the semiclassical approx- 
imation is not exact. In weak coupling, it leads to a scat- 
tering rate oc T, rather than cx T"^ j Ep, and does not give 
a mass renormalization; in strong coupling, while it gives 
correct Mott insulator behavior, i.e., divergence of the 
on-site self-energy cx l/(ici;„), the self-energy at small lo 
remains finite at T — > even in a paramagnetic metallic 
phase as long as the effective potential V{i^) has more 
than one degenerate minimum, implying incorrect non- 
Fermi liquid behavior. This originates from the neglect 
of the quantum fluctuation of HS fields. Including the 
quantum fluctuation to recover the correct Fermi-liquid 
behavior is not easyj^SiiMi Although, it is not exact, it 
will be shown below that the semiclassical method gives 
good estimates of the important physical quantities. 



FIG. 4; Example of the effective potential V (upper panel) 
and charge field ^ (lower panel) as functions of ip in the metal- 
lic Hubbard model on an infinite-dimensional FCC lattice 
[free DOS is given in Eq. @] with [/ = 4 and n = 0.5. Solid 
lines: T = 0.08 (about 14 % above Tc ~ 0.073); broken lines: 
T = 0.06 (about 18 % below Tc). The Hartree contribution 
(= nU) has been subtracted from ^. At both temperatures, 
potential is seen to be highly deviated from simple parabola. 
(/3-dependence of ^ indicates strong coupling between the two 
fields. 



There remains only one variable Lp, which is purely real. 
Numerical integrals can be performed without difficulty. 
Figure ^ shows the example of V{(p) and calculated 
by "full" semiclassical approximation for a non-integer 
filling. It is seen again that V is highly non-parabolic. 
Furthermore, depends on (p indicating the strong cou- 
pling between spin- and charge-fields. 



Now the approximate partition function Z, 



approx 



is ob- 
tained, one can obtain physical quantity from the func- 
tional derivative form following the DMFT procedure. 
The impurity Green function is 



^^^^ approx 

a = f 



1 



(18) 



where (...) stands for jdipe'x:p{—(3V(}p)}.../Zapprox- 
Thus, Eqs. (Unj-lO with the DMFT self-consistency 
equation construct the present semiclassical approxima- 
tion. 



III. WEAK COUPLING GAUSSIAN 
APPROXIMATION: COMPARISON OF 
CLASSICAL AND EXACT RESULTS 



This section uses the Gaussian approximation to the 
paramagnetic phase of the functional integral to gain in- 
sight into the limits of validity of the classical approxi- 
mation. The advantage of the Gaussian approximation 
is that any desired quantity can be computed straight- 
forwardly, permitting a detailed comparison of the re- 
sults of the classical approximation to the fully quantal 
calculation. We consider the mean square amplitude of 
the spin Hubbard-Stratonovich field, (t/?^), and also the 
Matsubara-axis electron self-energy, 5](ia;„). 

In the Gaussian approximation one determines the 

values Cp — UTTv^^,cr^z [ao + \{'f<^z + and x — 

iUTTi^^^^^ [ao -|- \ ((p<yz + is)] ^ which extremize the ar- 
gument of the exponential, and then expand the argu- 
ment of the exponential to second order in the deviations 
of the fields from their extremal values: 



Z Z f I?[cpa;]exp 



_<?(2) 



0(2 



(19) 
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where Z is contributions from the saddle points, and 



X;;. = ^'5(n-r2) 



i^Tr 



1 1 -1 

Go(t)== ao(T)l + -((^cr^ +xl)J(r)J , 



(20) 

,(21) 
(22) 



with dip = (7^ and 0:2: = 1, these are 2x2 matrices acting 
on spin space. In the mean-field approximation, a critical 
U exists; for U < Uc, ip — 0, and for U > Uc, ip ^ 0. 
(This transition is removed by fluctuations.) For U ^ Uc, 
spin fluctuations are very soft. 

We focus here on the ip integral, and we consider the 
small-?7 regime in which = 0, and in the paramagnetic 
phase, the x and (p integrals decouple. The mean square 
value of the fluctuations of ip, (t/?^) (obtained by perform- 
ing the Gaussian integral over all Matsubara components 
of ^p) and the classical approximation {v"^) ^lass (obtained 
by performing the integral only over zero Matsubara com- 
ponents of ip) are 



+ Uxo{m) ' 



-T- 



1 + C/Xo(0)' 

where xo is an irreducible susceptibility given by 



(23) 
(24) 

(25) 



with vi being a bosonic Matsubara-frequency. Finally, 
the leading perturbative contribution to the electron self- 
energy E and its classical approximation Y^ciass are given 
by 

T.{iuj„) = ^T'^xim) ao^iiujn -m), (26) 
EciassiiuJn) ^ ^TxiO) aQ^{iuJn), (27) 



with the susceptibility 



l + Uxo{m)' 



(28) 



We now present results obtained by applying the 
Gaussian approximation to the paramagnetic phase of 
the infinite-dimensional FCC lattice [DOS is given by 
Eq. In this section, energy is scaled by the vari- 

ance of DOS (==1). The lower band edge is at emin — 
— I/V2 « —0.71 and the square root divergence of the 
density of states means that for small fillings the density 
of states is very close to the bottom of the band. In the 
non-interacting limit the chemical potential is /i„=o.5 ~ 




FIG. 5: Ratio of the classical approximation to the full (Gaus- 
sian) mean square fluctuation of spin Hubbard-Stratonovich 
field i?x = {^^) class I i^'^) ^'^^ infinite-dimensional FCC lat- 
tice at densities n — 0.5 (upper panel) and n — 0.75 (lower 
panel). The curves cross because of the temperature depen- 
dence of the critical interaction strength (evident from the U 
at which the ratio approaches unity). 



—0.6335 {fi — Emin ~ 0.078) and ^,1=0.75 = —0.5352 
(m ^ Smm ~ .175). The square root divergences mean 
that the critical interaction strengths beyond which the 
mean field solution becomes unstable are themselves tem- 
perature dependent, and are small. 

Figure |S1 gives the comparison of the full and Gaus- 
sian fluctuation approximation to the mean square fluc- 
tuation of the spin Hubbard-Stratonovich field for the 
infinite-dimensional FCC lattice model with n = 0.5 and 
n = 0.75. One sees that the classical approximation cap- 
tures most of the fluctuations either for U near the crit- 
ical value or for temperatures of the order of the band- 
width. The panels of Fig. |H| show the full self-energy 
and the classical approximation for different interaction 
values, for n = 0.75. For small U, the semiclassical ap- 
proximation fails, but as U approaches the critical value 
(here approximately 0.84) the self-energy becomes well 
represented by its classical approximation 

The panels of Fig. [7| show the dependence of the self- 
energy on temperature, at a fixed moderate interaction 
strength (about 3/4 of the critical value). We see that 
at all temperatures studied, the classical approximation 
provides a reasonable estimate of the low frequency self- 
energy. At temperatures less than about half of the 
band width, the classical approximation grossly under- 
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FIG. 6: Frequency dependence of full (Gaussian) self-energy 
(filled circles) and classical approximation (open circles), for 
infinite-dimensional FCC lattice with n = 0.75, T = 0.04 and 
U values indicated. 



0.00 



30000000000000C 



0.00 




0.0 



-0.5 



.§"" -1.0 
w 

S -1.5 
-2.0 



- 






Ult = 6, Tit = 1/2 




— • — Semiclassical 




-O-QMC 



t -2 
-3 



-60 



10 



15 



20 



10 



15 



20 



10 



15 



20 



ft) It 



25 









Ult=S,Tlt=\l2 




— • — Semiclassical 




— O— QMC 







25 









;;///=20,m= 1/3 




— • — Semiclassical 




— 0— QMC 



25 



FIG. 8: Self-energies for a square-lattice half-filled Hubbard 
model computed from the semiclassical approximation (filled 
circles) and QMC (open circles) for interactions and temper- 
atures indicated. 



FIG. 7: Frequency dependence of imaginary part of full 
(Gaussian) self-energy (filled circles) and classical approxi- 
mation (open circles) for n = 0.75, U = 0.5 and T values 
indicated. 



estimates the high frequency part of the self-energy, but 
by T = 0.08 it is within about 30 % of the exact value, 
and for higher temperatures or for interactions close to 
the critical value the approximation is quite good. 



IV. COMPARISON OF SEMICLASSICAL 

APPROXIMATION TO QMC: 
TWO-DIMENSIONAL SQUARE LATTICE 

The Hubbard model on a bipartite lattice is known to 
exhibit an antiferromagnetic Neel ordering at half-filling 
and finite U. In this section, we apply the semiclassi- 
cal approximation to the Hubbard model with a square 
lattice [non-interacting electron dispersion is given by 
Eq. (Q], and investigate the self-energy, density of states 



and magnetic transition temperature. Note that the fi- 
nite T/v in two-dimensional square lattice is an artifact 
arising from the neglect of the low-lying spin-wave exci- 
tation which DMFT method can not capture. 

These results are compared with QMC performed 
on computing facilities at Universitat Bonn using the 
Harsch-Fye algorithm.* Typically, 48 time slices were 
used, and 2-20x10^ MC configurations were recorded. 
Occasional runs with more time slices or more MC con- 
figurations were made to verify error. Between 5 and 
(for the largest U) 40 iterations were needed for conver- 
gence of the DMFT loops. The Neel temperatures were 
determined from computations of the staggered magneti- 
zation. For U ^ 8t and T in or near the magnetic phase 
boundary, global update techniques (to be described else- 
where) were needed to ensure equilibration of the MC cal- 
culation. Error bars are smaller than the symbols shown. 

As noted in the previous section, it is expected that 
the semiclassical approximation becomes accurate in the 
strong-coupling regime and high temperature. In order to 
confirm this expectation beyond the harmonic (Gaussian 
fluctuation) approximation, we first calculate the elec- 
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FIG. 9: Density of states of square-lattice half-filled Hubbard 
model computed using semiclassical method (heavy lines) and 
QMC (light lines) for paramagnetic [panels (a) and (c)] and 
antiferromagnetic [panels (b) and (d), T ~ Q.KTm] phases at 
interaction and temperatures indicated. In the antiferromag- 
netic phase, majority spin density of states is shown by solid 
line, minority spin by broken line. 



tronic self-energies in the paramagnetic state and com- 
pare them with the QMC results. Recall that for this 
model the full band- width is 8i. In Fig.|Sl shown are the 
self-energies at J7/t = 6 and T/t — 1/2 (upper panel), 
U/t = 8 and T/t =1/2 (middle panel), and U/t = 20 
and T/t — 1/3 (lower panel) computed by the semiclas- 
sical approximation and QMC as functions of imaginary 
frequency. With increase of U , the self-energy at low 
frequency increases, and for U > Uc, diverges indicat- 
ing the Mott metal-insulator transition. This behavior is 
well reproduced by the semiclassical approximation, and 
the agreement with the QMC is remarkable, particularly 
in the strong coupling regime. 

Figure compares the semiclassical results for a real- 
frequency quantity, the density of states, to results ob- 
tained by maximum entropy analytical continuation of 
QMC data. This is a rather stringent rest of the method, 
and agreement is seen to be reasonably good. In the weak 
coupling, paramagnetic phase [Fig. 0(a)], the semiclas- 
sical approximation underestimates the lo = Q peak (be- 
cause a; = scattering rate is overestimated) and under- 
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FIG. 10: Neel temperature of square-lattice half-filled Hub- 
bard model as functions of U computed from semiclassical ap- 
proximation (filled symbols) and QMC (open symbols), HF 
(hght solid line) and large U limit {Tn = 4t'^/U) (light broken 
line) . 



estimates the weight in the wings (because the high fre- 
quency scattering rate is underestimated, see the upper 
panel of Fig. IS)). Below Tn, the magnetization is slightly 
overestimated by the semiclassical approximation (mag- 
netization (to) ~ 0.54 for the semiclassical approximation 
and (to) - 0.4 for QMC at U/t = 4 and T/t = 0.2) as can 
be seen from the slightly larger gap. DOS is sharper be- 
cause the self-energy at high-frequency is underestimated 
as in the paramagnetic state. In the strong coupling 
regime, agreement between the semiclassical approxima- 
tion and QMC is quite good as can be expected from the 
self-energy (see the lower panel of Fig.lSJ. 

Figure IIUI summarizes results for the magnetic phase 
diagram obtained from the semiclassical approximation, 
QMC, the HF approximation and a strong-coupling ex- 
pansion. The semiclassical approximation is seen to give 
remarkably good results over the whole phase diagram. 
In the weak-coupling limit, the semiclassical approxima- 
tion gives correct behavior: Tn asymptotes to the re- 
sult of HF in the limit of U 0. This is natural be- 
cause the two approximations becomes identical in the 
[/ — > 0, T ^ limit. At finite but small U region, the 
reduction of Tn from HF is found to be quite large, be- 
cause the finite self-energy reduces nesting effect substan- 
tially. The present approximation gives correct behavior 
in the strong-couphng regime; Tn asymptotes to At^ /U, 
the mean-field result of the strong coupling expansion, 
and, in particular, is almost identical to the QMC re- 
sults. It should be noted that at large U and low T, 
the QMC requires a very large amount of Monte-Carlo 
sampling to reach equilibrium, whereas the semiclassical 
method is numerically very cheap. Even for U/t — 16 
and T/t ~ 0.2 (antiferromagnetic phase), it takes less 
than 5 mins. on a commercial PC to compute one point 
in the present method, while it takes ~ 6 hrs. by QMC 
on a similar computer for the same parameters. 
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V. COMPARISON OF SEMICLASSICAL 
APPROXIMATION TO QMC: FCC LATTICE 

In this section, we apply the semiclassical approxima- 
tion to the single-band Hubbard model on the FCC lat- 
tice in infinite- and three-dimensions, and compare the 
results to QMC data of Ulmkc,^! and to the "two-site" 
approximation of Potthoff.'^'^ We focus on the filling de- 
pendence of the magnetism in these models, in which 
the charge fluctuation plays an important role. At non- 
integer flUing, these models order ferromagnetically be- 
cause of the large DOS near the bottom of the band, 
in contrast to the antiferromagnetic ordering found in 
half-filling square lattice. In addition, near half-filling, 
the three-dimensional model is found to show an addi- 
tional competing order: a "layer-type antiferromagnetic" 
state. We examine whether the semiclassical approxima- 
tion capture this behavior. In this section, we use the 
variance of DOS v as an unit of energy. 

First, we apply the semiclassical approximation to 
the infinite-dimensional FCC lattice with DOS given in 
Eq. (Q, V — 1. The density of states above the Curie tem- 
perature is essentially the same as seen as the solid line in 
Fig. 121 and consists of two structures: one peak just below 
the Fermi level a; = and shoulder at w ~ nU . These cor- 
respond to the lower- and upper- Hubbard bands. (Note 
that the result shown in Fig. O is computed at T < Tc, 
but in the paramagnetic phase, by suppressing the mag- 
netic transition.) The separation between the lower- and 
upper-Hubbard band is somewhat reduced from the bare 
value of U (by about a factor of mean occupation num- 
ber n), possibly because the neglect of the fluctuation 
of a charge-fleld. The two structures continuously evolve 
into the majority-spin and minority-spin bands below Tc 
as shown in Fig. ITTI 

Figure [TTl also presents the DOS from QMC. (Temper- 
ature and density are chosen so that the magnetization 
in the two calculations are similar.) In QMC results, we 
observe a sharp peak below Fermi level for the majority 
spin band. For the minority spin band, there appears 
a sharp peak near the Fermi level and a broad hump 
at w ^ 3. The former originates from the fact that 
the system is not fully polarized, and the latter corre- 
sponds to the upper-Hubbard band. As noted above, the 
upper-Hubbard band (broad hump around w = 3 in the 
minority-spin band) appears slightly low in energy in the 
semiclassical approximation. Except for this discrepancy, 
the overall structure of DOS is well reproduced by the 
semiclassical approximation. In particular, the occupied 
bands below lo ~ Q show reasonable agreement. 

In Fig. 1121 ferromagnetic Curie temperatures Tq for 
different values of U are shown as functions of electron 
density. Filled symbols are the results obtained using 
the semiclassical approximation. For comparison, results 
of QMG^ and HE are shown as open symbols and light 
lines, respectively (Note that the HE Tc shown are re- 
duced by a factor of 10) . The semiclassical approximation 
overestimates Tc than QMC by a factor of ~ 50 % over a 
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FIG. 11; Comparison of the density of states of infinite- 
dimensional FCC Hubbard model with f/ = 4 and temper- 
atures chosen so that (m) = 0.4 in each case. Heavy lines: 
results of the semiclassical approximation at n = 0.6 and 
T = 0.07 (~ 13 % below Tc ~ 0.08). Light lines: results of 
QMC at n = 0.58 and T = 0.04 (~ 20 % below Tc ~ 0.05) 
taken from Ref.l41i Solid and broken lines are for the majority 
and minority spin. 

wide range of density. (Of course, the difference near the 
critical density of QMC is much larger.) It is seen that 
HE approximation is very poor in all parameter regimes, 
highly overestimating the ferromagnetic Curie tempera- 
ture. The critical density nc is also overestimated, at 
Tie ~ 1.2 at [/ = 2 and - 1.5 at C/ = 4. The higher Tc 
found in the semiclassical approximation may be due to 
the neglect of the quantum fluctuation of the HS fields. 
Reduction of the transition temperature due to the quan- 
tum fluctuation can be seen in the context of electron- 
phonon coupling in Ref. However, good agreement 
between the semiclassical approximation and QMC indi- 
cates that the thermal part of the fluctuation dominates 
the electron self-energy in the wide region and, thus, the 
magnetic transition. Critical densities ric where the fer- 
romagnctism disappears at T = are found to be slightly 
higher in the semiclassical approximation than in QMC, 
nc 0.88 for [/ = 2 and - 0.97 for = 4, while 
QMC gives - 0.7 for [/ = 2 and - 0.88 for U = A. 
However, in the light of the simplification of the present 
approximation, the agreement with QMC within 20 % 
error is remarkable. 

We applied the two-site DMETiS^ to the same model 
to investigate the magnetic phase diagram. The critical 
densities ric obtained by two-site DMET are found to be 
very similar to those by the semiclassical approximation. 
However, Curie temperatures are found to be overesti- 
mated by a factor of ^ 6 compared with QMC, and by 
a factor of ~ 3-4 compared with the semiclassical ap- 
proximation. In the two-site DMET, Tc at n = 0.5 is 
found to be - 0.17 for U = 2, and ~ 0.26 for [/ 4. 
The overestimate arises because, in the two-site DMET, 
the lAM is composed of only two sites, one is correlated 
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FIG. 12: Ferromagnetic Curie temperature of single-band 
Hubbard model on an infinite-dimensional FCC lattice as 
function of electron density for U = 2 (squares) and 4 (cir- 
cles). Filled symbols: semiclassical approximation; open sym- 
bols: QMC (Ref.0). Curie temperatures (xlO~^) computed 
by HF approximation to the same model with U — 2(4) are 
shown as a light solid (broken) line. Two-site DMFT results 
for the transition temperatures are not shown, but represen- 
tative two-site results for n = 0.5 are Tc ~ 0.17(0.26) for 
U = 2(4). Phase boundary at T = for [/ = 2(4) computed 
by the two-site DMFT is shown as a filled (open) triangle. 



(impurity) and one is non-correlated (bath). The energy 
difference between these sites is generically large, so the 
(spin) entropy is underestimated. The two-site DMFT 
also yields reentrant behavior near the critical concentra- 
tion, which is not observed in the semiclassical approxi- 
mation and QMC. Adding more sites in the "bath" part 
of 1AM would remedy these behaviors, but at drastically 
increased computational expense. 

Finally, we investigate the magnetic instability in the 
more realistic three-dimensional FCC lattice whose free 
band-dispersion is given in Eq. |(5J- Following Ref. 0, 
we take t' = t/A and choose the variance of the DOS 
V — \/ 12t^ + 6t'^ as the energy unit, thus t « 0.284 and 
t' « 0.071. With this parameter set, the bottom of the 
band is given by £mm = — 4i + 2t' « —0.994. There is 
no divergence in the free DOS in contrast to the infinite- 
dimensional case. 

The magnetic phase diagram of the three-dimensional 
FCC lattice computed by the semiclassical approxima- 
tion is shown in Fig.^1 For comparison, QMG^ and HF 
results are also plotted. It is seen that the semiclassical- 
approximation results for Tc are about twice higher than 
QMC, while HF results are about 20 times higher than 
QMC. Similar to the infinite-dimensional case, higher 
critical density Uc is overestimated by ~ 20 %. QMC 
calculations are argued to indicate that there exists a 
lower critical density n'^ ~ 0.15 below which the fer- 
romagnetism disappears^! This behavior is ascribed to 
the absence of the divergence at the bottom of the bare 
DOS. In the semiclassical approximation, Tc seems to be- 
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FIG. 13: Upper panel: Phase diagram of single-band Hub- 
bard model on a three-dimensional FCC lattice as a function 
of electron density and temperature for interaction ?7 = 6 and 
variance of the density of states v = 1. Squares and circles are 
Curie temperature and Neel temperature for layer-type anti- 
ferromagnetic state, respectively. Filled symbols: results of 
the semiclassica l ap proximation: open symbols: QMC results 
taken from Ref.l4ll Curie temperature (xlO"'^) computed by 
HF approximation to the same model is shown as a light solid 
line. Lower panel: Expansion of the phase diagram, showing 
the region near n = 1 enclosed by a broken line in the up- 
per panel. The semiclassical results for Tc and Tjv indicate 
the instability of ferromagnetic and antiferromagnetic states, 
respectively, towards the paramagnetic state. 



come zero proportional to the electron density similarly 
to the HF approximation. In the latter, clearly n'^ — 0. 
This might be because the semiclassical approximation 
reduces to the HF at weak coupling (in this case small 
density) and low temperature limit. However, it is very 
difficult to judge if n'^ = or not in the present accuracy 
for the semiclassical approximation. 

Differently from the infinite-dimensional FCC lattice, 
QMC calculation found another magnetic state stabilized 
near n = 1: a layer- type antiferromagnetic state with a 
magnetic vector q = (tt, 0, 0).'^^ The semiclassical approx- 
imation is successful in finding the layer-type antiferro- 
magnetic state near n — 1. Computed Neel temperature 
at n = 1 is Tat ~ 0.048 which agrees with the QMC 
result within the statistical error of QMC. Better agree- 
ment in the Neel temperature than in the Curie tem- 
perature may be attributed to the large U used in this 
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calculation and to the suppression of the charge fluctua- 
tions near half-fiUing. According to the two-site DMFT, 
critical Uc to the Mott transition at half filling is esti- 
mated to be Uc = 6v (correct value is expected to be 
slightly smaller then this).^'^ With the parameter used, 
U = 6, the system is in a Mott insulating state at n = 1, 
and the charge fluctuation is suppressed. Therefore, the 
thermal fluctuation of spin field dominates the magnetic 
transition. In contrast to QMC, ferromagnetic and anti- 
ferromagnetic states contact with each other in the semi- 
classical approximation. (In the semiclassical result, Tq 
and T^r indicate the instability of the ferromagnetic and 
antiferromagnetic states, respectively, to the paramag- 
netic state.) In the light of the overestimation of upper 
critical density Uc for the ferromagnetism, this failure is 
also supposed to be from the neglect of the quantum fluc- 
tuation of HS fleld. This point remains to be addressed 
in the future work. 

In the metallic region, one needs to flx the chemical 
potential according to the density- this must be done at 
each interaction strength and temperature, which is time 
consuming. Even in this case, the semiclassical scheme 
is found to be computationally very cheap. Typical CPU 
time is less than 5 mins. a.t U — 2 and T = 0.05 using a 
commercial PC. 

Summarizing this section, we investigated the mag- 
netic behavior of the Hubbard model on infinite- and 
three-dimensional FCC lattices using the semiclassical 
approximation. The magnetic phase diagram com- 
puted by the semiclassical approximation show reason- 
able agreement with QMC. The ferromagnetic Curie tem- 
perature in metallic region is found to be within a factor 
of '--^ 2 compared with QMC, and the antiferromagnetic 
Neel temperature near n = 1 is found to agree with QMC 
within statistical error. Better agreement between the 
semiclassical approximation and QMC in the Neel tem- 
perature near integer-flUing is supposed to be from the 
fact that the charge fluctuation is suppressed by the cor- 
relation, so that thermal spin fluctuations dominate the 
transition as in the half-filled square-lattice case. Poorer 
agreement in the Curie temperature is expected to be 
from the neglect of quantum fluctuation or charge fluc- 
tuation. 



VI. APPLICATION OF SEMICLASSICAL 
APPROXIMATION TO DCA AND 
FICTIVE-IMPURITY METHOD 

Our semiclassical approximation is easily combined 
with cluster schemes, such as dynamical cluster 
approximation^^ (DCA) and real-site cluster extension 
of DMFTi^S^Siii As a simplest example, we use the two- 
site DCA and fictive-impuritjiS^ (FI) methods to study 
the Hubbard model on a square lattice. 

We consider half-filling, thus charge fields are absorbed 
into the chemical potential shift. The Weiss field becomes 
2x2 matrix a, and in the paramagnetic phase, this has 
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FIG. 14: Nearest-neighbor spin correlations —{o'izcr2z) of a 
square-lattice Hubbard model {U/t = 20) as functions of 
T computed using two-site DCA (open circles) and the FI 
method (open squares), compared to QMC (filled symbols) 
and high-temperature series expansion(light line). 



a form 



(29) 



with ao and oi being the on-site and NN Weiss fields, 
respectively. 1 and are the 2x2 matrices acting on 
orbital (impurity site) space. The HS transformation is 
performed at each impurity sites i = 1,2 with spin field 
(fi. The effective potential for HS field becomes 



^^(¥^1,^2) = ^(V'l ^ 
-TTrln 



a 



^2 



(30) 



Here, cr^ is acting on spin space, and Tr is taken over 
Matsubara-frequency, orbital and spin indices. Then, the 
partition function is given as 



dipid(p2 ex.p{-f3V{ipi,ip2)}- 



(31) 



Finally, on-site and NN-site Green functions, Gq and Gi, 
respectively, are obtained via 



^ 1 ^ In Zfipp^'ox 

2 Sao ' 
1 (5 In Zfippj-ox 

^'^2 6^, ■ 



(32) 
(33) 



The self-consistency equations for DCA are closed follow- 
ing the scheme presented in Ref. The self-consistency 
equations for FI method and the general formalism for 
the cluster extension of DMFT are given in Ref. IbM 

Nearest-neighbor spin correlations —{aiz<J2z) com- 
puted by DCA and FI with the semiclassical approxima- 
tion and QMC as functions of temperature for U = 20t 
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are shown in Fig. 1141 Similar spin correlation obtained 
by QMC in both DCA and FI supports the applicability 
of the semiclassical approximation to the cluster DMFT. 
For comparison, the same quantity computed using a 
high-temperature series expansion (HTS) for the = 1/2 
NN Heisenberg model with J = /U are also shown. It 
is seen that spin correlation obtained from the DCA is 
much larger than the HTS result. We believe the origin 
of the discrepancy between DCA and HTS is that DA 
is equivalent to imposing a periodic boundary condition 
in all the directions in a real-space cluster is adopted. 
Thus, in the two-site DCA, one has a model in which 
two sites are connected via z bonds [connectivity z — A 
in a square lattice] and the spin correlation is thus over- 
estimated. Interestingly, FI method gives almost identi- 
cal curves to HTS. (slight deviation can be seen below 
T ^ QAt.) A detailed study of multi-site cluster models 
using semiclassical methods and QMC will be presented 
elsewhere 4^ 



VII. EQUILIBRATION AND PARTITION 
PHASE SPACE 

In this section, we point out an additional advantage of 
the semiclassical approximation. A key issue in Monte- 
Carlo simulations is equilibration, which is particularly 
difficult in systems in which the phase space is partitioned 
into several nearly equivalent minima, separated by large 
barriers. Local update techniques require extremely long 
runs to climb over barriers, while global update tech- 
niques are expensive and sometimes inconvenient to im- 
plement and bring their own convergence issues. The 
partitioned phase space phenomenon occurs frequently 
in strongly-correlated models, and represents a signifi- 
cant obstacle to practical computations. The semiclas- 
sical approximation, by contrast, is inexpensive enough 
that the entire (semiclassical) phase space can be sam- 
pled. 

For example. Fig. El shows the distribution func- 
tion of spin HS fields for the two-impurity DCA for 
the square-lattice half-filled Hubbard model defined by 
P{ipi,Lp2) = exp{~(3V{ipi,ip2)}/Z approx computed at 
temperature T/t = 0.5 > Tn (see Fig. The par- 

titioning of phase space is evident. It is seen that 
the distribution has twofold symmetry, not fourfold, in 
ipi — ip2 plane; ipi > ip2 and ipi < ip2 are not equiv- 
alent. Larger peaks at (1^91,(^2) — (±20t,^20t) than 
((/3i,(^2) — (±20t,±20t) indicate that the antiferromag- 
netic correlations prevail far above the Neel temperature, 
giving the result shown in Fig. E| Reliable estimates of 
— (o'i2CT2z) requires sampling of all four extreme, with cor- 
rect weights, which is very difficult to achieve in QMC at 
low T (this is why we do not present data at T < 0.5i, 
and statistical errors are evident even at T = 0.6f). How- 
ever, the integrals involved in the semiclassical method 
may be performed with no difficulty. 




FIG. 15: Distribution of the spin fields ip\ and <^2, P['fi, '^2), 
computed by two-impurity DCA for the Hubbard model 
on a square lattice. P((pi,Lp2) is defined by P{(fii,ip2) = 
exp{~f3V{ipi,(p2)}/Zapprox. Parameters are U/t — 20 and 
T/t = 0.5. Larger peaks at {•^1,1^2) ^ {±20t,=f20t) than 
((^1,(^2) — (±20f, ±20t) indicate antiferromagnetic correla- 
tion. 



VIII. SUMMARY AND DISCUSSION 

In this paper, we investigated the semiclassical approx- 
imation to the continuous Hubbard- Stratonovich trans- 
formation as an impurity solver of DMFT method for the 
correlated-electron models. The Hubbard-Stratonovich 
transformation introduces two auxiliary fields, coupling 
to the electron spin and charge density respectively. The 
semiclassical approximation consists of retaining only the 
classical (zero-Matsubara-frequency) component in the 
functional integral over the spin field, and, for each value 
of the spin field, using a steepest descents approximation 
to approximate the integral over the charge field by the 
value which extremizes the action at the given value of 
the spin field [see Eq. H15|) ]. This treatment of the charge 
field was found to be essential for achieving reasonable 
results: see Fig. |31 and the associated discussion. The 
semiclassical approximation captures the thermal fluc- 
tuation of spin field efficiently beyond harmonic (Gaus- 
sian fluctuation) approximation and is applicable to both 
the metallic and insulating region. Estimates obtained 
using the Gaussian fluctuation approximation indicate 
that the semiclassical approximation gives reasonable re- 
sults at larger U (when spin fluctuations are soft) or at 
high temperature. (Sec. Ill) We applied the semiclassical 
approximation to the single-band Hubbard model on a 
two-dimensional square lattice (half-filling) (Sec. IV) and 
infinite- and three-dimensional FCC lattices (finite dop- 
ing) (Sec. V). The semiclassical approximation is found 
to give reasonable results, with accuracy improving for 
stronger couplings and in situations where charge fluctu- 
ations are suppressed. A particularly attractive finding 
is that the procedure finds multiple phases when these 
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exist. (Fig. I13|l The key point appears to be that, at 
stronger correlations, the physically important fluctua- 
tions become very soft (justifying a semiclassical approx- 
imation) but very anharmonic (requiring an integral over 
all field configurations). Due to the neglect of quantum 
effects, the semiclassical approximation fails to reproduce 
correctly the quasiparticle resonance, which will certainly 
be important at weak coupling or at low temperature. 
One possible improvement of this failure would be us- 
ing a interpolative method^^ to reproduce the low-energy 
quasiparticle in a paramagnetic metallic region. For the 
half-filled square lattice, we showed the semiclassical ap- 
proximation gives reasonable behavior of the self-energy 
which is comparable to QMC in a wide range of param- 
eters. The density of states and Neel temperature com- 
puted by the semiclassical approximation are found to be 
in very good agreement with QMC. In the metallic FCC 
lattice, the DOS has a two-peak structure correspond- 
ing to the lower- and upper-Hubbard bands in a para- 
magnetic phase. These structures evolve into majority- 
and minority-spin bands as the temperature is decreased 
through the ferromagnetic transition. Comparison of the 
DOS computed by the semiclassical approximation and 
by QMC shows reasonable agreement, especially as con- 
cerns the occupied state. Ferromagnetic Curie tempera- 
tures computed by the present method are found to be in 
the range of factor of 2 compared with the QMC. T = 
phase boundaries are found to be within the error of 20 %. 
The semiclassical approximation for the Curie tempera- 
ture shows better agreement than two-site DMFT, and 
T — phase boundaries are found to be in the same 
range of two-site DMFT. As for the three-dimensional 
FCC lattice, the semiclassical approximation is able to 
detect antiferromagnetic state near the integer filling con- 
sistent with QMC. Neel temperature in this case agrees 
with QMC within the statistic error. In Sec. VI, we 
apply the semiclassical approximation to two-impurity 
DCA and real-space cluster DMFT (fictive-impurity) for 
the square-lattice Hubbard model at half-filling, and in- 
vestigate how spatial correlation is taken into account in 
the present approximation. We confirmed that the short- 
range (nearest-neighbor) spin correlation prevails above 
the Neel temperature. Comparison of nearest-neighbor 
spin correlation between the semiclassical approximation 
and QMC shows good agreement in ranges where the 
QMC calculations are well converged. In Sec. VII, equi- 
libration problem in QMC associated with the partition 
of phase is discussed. This problem does not occur in the 
semiclassical approximation. 

In conclusion, we mention some open scientific ques- 
tions which can be addressed by the semiclassical 
method, and mention two areas in which further inves- 
tigation and improvement would be desirable. A key 
opportunity involves systems with several sites and/or 
several orbitals. In such cases, QMC suffers from "sign" 
problems, and both QMC and ED method scale poorly 
with system size. The semiclassical method does not 
suffer from "sign" problem; while there appeared sim- 



ilar problem associated with the imaginary coefficient 
for the charge-field, this is resolved by applying the 
saddle-point approximation to the charge-field. When 
it is applied to systems with > 1 orbitals (on one 
or more sites), computational time for the semiclassical 
approximation scales of with P^p and A^ being a to- 
tal number of descretized HS field and total number of 
sites/orbitals, respectively. By contrast in Monte-Carlo 
methods, the scaling is [{(3U)'^Nmc]^ (A^a/c is a number 
of Monte-Carlo samphng, typically Nmc ~ 10^), and by 
4^Nii+N,^tH) for ED {Ntath is a number of "bath" sites 
per impurity orbital, and Ntath should be larger than 
2-3 from the overestimation of Tc by two-site DMFT). 
As can be seen in the Fig. 1151 dominant contribution is 
known to be from \ip\ U , and if need it would be possi- 
ble to simplify the A^-dimensional integral over HS fields. 
Thus, the semiclassical approximation appears to be an 
attractive option for studying larger clusters than QMC. 

As for the multisite problem, we have only applied the 
semiclassical approximation to the single-orbital Hub- 
bard model at half-filling with DCA and FI method. In 
this case, the charge-field is absorbed into the chemical 
potential. Away from the half-filling, one needs to fix 
the charge-field at each configuration of the spin-fields 
according to Eq. IjlSfl generalized to the multisite situa- 
tion. There is no proof that Eq. H15|l has unique solution 
in this situation although we have so far not encountered 
problems. This issue needs further investigation. 

As one of the applications of the semiclassical ap- 
proximation, investigation of the systems with degen- 
erate orbitals like transition-metals oxides with Cg or 
t2g electrons would be interestingiSi^ Theoretical stud- 
ies of the ferromagnctism in these situation have been 
presented (Refs. 33 34-3 j) . but these works dealt with 
spin-fluctuations in the metallic state, fluctuation and 
ordering associated with the quadrupole moment for Cg 
and/or t2g orbital was neglected. Hubbard-Stratonivich 
transformation including spin and quadrupole moment 
for two-band Hubbard interaction for Cg orbital has been 
introduced in Ref. l3^ However, pin and orbital transi- 
tion at finite temperature using the multi-band Hubbard 
model with full symmetry—^ remains to be investigated. 

Another promising application of the semiclassical ap- 
proximation would be the spatially inhomogeneous sys- 
tems, where the computational expense associated with 
other aspects of the problem renders an inexpensive im- 
purity solver essential. Recently, two of the authors (S.O. 
and A.J.M.) apphed two-site DMFT to the spatially- 
inhomogeneous heterostructure problem and investigated 
the evolution of the electronic state (quasiparticle band 
and upper- and lower-Hubbard bands) as a function of 
distance from the interfacei^ Applying the semiclassical 
approximation to such systems to investigate the possi- 
ble spin (and orbital) orderings^ is an interesting and 
urgent task. 
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